The high-energy collision of two black holes 
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We study the head-on collision of two highly boosted equal mass, nonrotating black holes. We 
determine the waveforms, radiated energies, and mode excitation in the center of mass frame for a 
variety of boosts. For the first time we are able to compare analytic calculations, black hole per- 
turbation theory, and strong field, nonlinear numerical calculations for this problem. Extrapolation 
of our results, which include velocities of up to 0.94c, indicate that in the ultra-relativistic regime 
about 14 ± 3% of the energy is converted into gravitational waves. This gives rise to a luminosity 
of order 10 _2 c 5 /G, the largest known so far in a black hole merger. 
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/. Introduction. An important and long-standing prob- 
lem in general relativity concerns the ultra-relativistic 
scattering of black holes (BHs). This is one of the most 
violent events one can conceive of in the theory. The lack 
of solutions has spurred much speculation about what 
may happen in this regime. For example, these events 
are a natural testing ground for the cosmic censorship 
conjecture: is there a class of initial conditions where 
they generically lead to the formation of a naked singu- 
larity, or do event horizons always form to cloth singular 
behavior in the geometry? 

Related questions concern the ultra-relativistic scat- 
tering of particles. If the center of mass (CM) energy is 
beyond the Planck scale, gravity is expected to dominate 
the interaction. Furthermore, since the kinetic energy 
dominates over the rest mass energy, the gravitational in- 
teraction should be rather insensitive to the structure of 
the particles, implying that the trans-Planckian scatter- 
ing of point particles should be well described by BH scat- 
tering [l| . This is of particular relevance for recent pro- 
posals to solve the hierarchy problem by adding "large" 
extra dimensions 0], or an extra dimension with a warp 
factor [3[ , thus producing an effective electroweak Planck 
scale. This offers the exciting possibility that BHs could 
be produced in particle colliders and ultra high-energy 
cosmic ray interactions with the atmosphere [l], 0]. A 
naive estimate of the cross section for Mp\ ~ 1 TeV pre- 
dicts that super- TeV particle colliders will produce BHs 
at a rate of a few per second, making the Large Hadron 
Collider (LHC) at CERN a potential black hole factory. 
An important element to search for BH production sig- 
natures is to understand the BH scattering process, and 
in particular the energy lost to gravitational radiation. 
Given that the beam commissioning to 7 TeV is sched- 
uled for late 2008, this is a timely research topic. Fur- 
ther interesting applications of high-speed BH collisions 
to high-energy physics have recently been suggested by 
the AdS/CFT correspondence conjecture Particu- 



larly intriguing is the possibility of using this duality to 
understand properties of the quark-gluon plasma formed 
in gold ion collisions at Brookhaven's Relativistic Heavy 
Ion Collider (RHIC) through a study of ultra-relativistic 
BH collisions in AdS 0. 

Early attempts to understand the ultra-relativistic BH 
scattering problem were based on work by Penrose [7[ in 
the 1970s. He modeled the spacetime metric as the union 
of two Aichelburg-Sexl waves [8j , describing the collision 
of two infinitely boosted Schwarzschild BHs, and found 
a closed trapped surface at the moment of collision, giv- 
ing an upper limit of roughly 29% of the initial energy of 
the spacetime radiated in gravitational waves. Beyond 
the collision event the solution is unknown. Given the 
extreme conditions of high-speed scattering it is unlikely 
that analytic solutions describing the full dynamics of the 
spacetime will be found, and therefore numerical meth- 
ods must be employed. Only recently have long-term 
stable numerical evolutions of black-hole binaries been 
achieved Q. The flurry of subsequent activity explor- 
ing the merger process has so far exclusively focused on 
rest-mass dominated scenarios (see [l(| for a review). 

In this Letter we report the first numerical solutions de- 
scribing the collision of two equal mass BHs in the regime 
where the initial energy of the system is dominated by 
the kinetic energy of the BHs. In Sec. II we describe the 
problem setup, including the numerical code and initial 
conditions. We also review some existing analytical ap- 
proximations to aspects of the problem, which will be 
important both to interpret the numerical results and to 
give some confidence in extrapolations of the results to 
infinite boost. In Sec. Ill we present the primary results, 
focusing on the gravitational waves emitted during the 
collision. Concluding remarks are given in Sec. IV. Un- 
less stated otherwise, we use geometrical units G = c = 1. 
//. Numerical Setup and Analysis Tools. The nu- 
merical simulations presented here have been performed 
with the Lean code, described in detail in [1JJ], where 
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head-on collisions of different classes of initial data were 
compared. Here we exclusively study evolutions of punc- 
ture initial data [13] describing two equal mass, non- 
spinning, boosted BHs colliding with zero impact param- 
eter in the CM frame. The initial coordinate separation 
between the punctures is set to ro, and the boosts are 
prescribed in the form of non- vanishing Bowen-York [l3j 
parameters ±P for the initial linear momentum of either 
BH. The Hamiltonian constraint is solved using Ansorg's 
spectral solver TwoPunctures The irreducible 

masses M- ml 2 of the BHs are estimated from their appar- 
ent horizon areas, calculated using Thornburg's apparent 
horizon finder AHFinderDirect This enables us 
to calculate the BH masses M1.2 from Christodoulou's 
[HI ] relation Mf 2 = M 2 ITl 2 + P 2 1 from which we define 
the Lorentz boost parameter 7 = Mi^/Mjai^ (cf. [III). 
From a numerical point of view, simulations with large 
values of 7 are challenging, partly because the Lorentz 
contraction decreases the smallest length scale that needs 
to be resolved. Thus mesh-refinement is essential, and 
here it is provided via the Carpet package [Lsl ]. 

We use the Newman-Penrose scalar ^4 to mea- 
sure gravitational radiation. At an extraction radius 
r from the center of the collision we decompose VP 4 
into multipole modes ipim of the spherical harmonics 
of spin-weight —2, -2^ m , according to ^(i, r, 6, 4>) — 

X^2l]L=-z -"2Yi m (8 ,4>)<4)i m {t,r). Due to the symme- 
tries of this problem, the only non-vanishing multipolcs 
all have even I, m = 0, and are purely real, corresponding 
to a single polarization state h + . The energy spectrum 
and luminosity of the radiation are given by 
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respectively, where a hat denotes Fourier transform. 

Our results are affected by three main sources of uncer- 
tainties: finite extraction radius, discretization and spu- 
rious initial radiation. We reduce the error arising from 
finite extraction radius by measuring the waveform com- 
ponents at several radii, and fitting them to an expression 
of the form ipi m (r, t) = ipj® (t) + ip^ (t) jr. The waveform 

"at infinity" ipj® (t) is the quantity reported throughout 
this work and used to calculate related quantities, such as 
the radiated energy. The uncertainty in this extrapolated 
value is estimated by performing a second fit including 
also a quadratic term ip^/r 2 , and taking the difference 
between the first- and second-order fits. The resulting 
uncertainty in the radiated energy is typically ~ 3 — 5 %. 

To estimate discretization errors we evolved the most 
challenging simulation with 7 rs 3 with resolutions h — 
M/174, M/209 and M/244, where M = M x + M 2 . We 
observe convergence slightly below second order in the 
total radiated energy, and use a conservative estimate of 



10% for the resulting error near 7 ~ 3, which drops to a 
few percent in the non-boosted case (cf. [Til]). 

Finally, the conformally flat puncture initial data is 
known to contain spurious gravitational radiation, which 
increases strongly with boost 7 (from a few times 10~ 5 
for BHs at rest, to about 8% of the total ADM mass of 
the system for 7 « 3). In order to extract physically 
meaningful information, one has to separate the spurious 
radiation from the radiation generated by the collision it- 
self. This is done by "waiting" for the spurious radiation 
to pass the last extraction radius, and then discarding 
the earlier part of the wave signal. For large boosts, the 
amount of time between the trailing edge of the spuri- 
ous radiation and the leading edge of the waves emitted 
during the collision is roughly ro/(47 2 ). Thus, the ini- 
tial separation required to cleanly extract the emitted 
signal increases rapidly with 7. Because large separa- 
tions require larger computational domains and longer 
run-times, the spurious radiation effectively limits our 
ability to study very large 7. With current resources, 
we were able to use initial separations of up to 66M for 
7 > 2, leading to an uncertainty in the total radiated 
energy which grows rapidly with boost, reaching a value 
of 5% for 7 w 3. By combining all errors, we estimate 
the total uncertainty in the radiated energy to be about 
15% for 7 ~ 3, about 10% near 7 = 2 and a few percent 
for simulations with small velocities. 

High-energy collisions are uncharted territory for nu- 
merical relativity. It is helpful, therefore, to have al- 
ternative methods for guidance and consistency checks. 
Besides Penrose's bound, we will make extensive use of 
extrapolations of Smarr's "zero-frequency limit" (ZFL) 
fl9| and of point particle (PP) calculations [2(J, where 
one considers a small object of mass m colliding with a 
massive BH of mass Mbh to linear order in to/A/bh- 
III. Results. We ran a series of simulations from 7 = 1 
to 7 ks 3, with initial separations as discussed in the pre- 
vious section. In all cases the collision results in a single 
BH plus gravitational radiation, i.e. there is no sign of 
any violation of cosmic censorship. The final BH is born 
highly distorted. We measure the distortion by taking 
the ratio C of the proper equatorial to polar circumfer- 
ences of the common apparent horizon (CAH). For the 
range of boosts studied here, the peak value is well fit- 
ted by the relation C pcak - 1.5 - O.5/7. Thus in the 
large-7 limit C peak ~ 1.5, in agreement with Penrose's 
result of C = 7r/2 for a CAH consisting of two flat disks. 
After birth, the BH settles down to a Schwarzschild so- 
lution, and the gravitational radiation can be described 
as a superposition of quasinormal modes (QNMs) of the 
resulting BH. 

In Fig. [T] we show the dominant component "020 of the 
waveform from collisions with 7 = 1.07, 1.3, 1.7, 3.0 
(corresponding to (3 = v/c ~ 0.36, 0.64, 0.82, 0.94, re- 
spectively). The origin of the (t — r) axis roughly corre- 
sponds to the instant of formation of a CAH. One can 
identify three main parts in the waveforms: a precur- 
sor, a main burst at the onset of the CAH formation 
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FIG. 1: Dominant multipolar component tp2o(t — r) for dif- 
ferent values of (3, as indicated in the inset. 




FIG. 3: Total radiated energy (including error bars) as a func- 
tion of /3, and best fit using the ZFL prediction. 



3 

T3 

CM 

LU 

"a 10" 











— P=0.94 
: (5=0-82 
: —- (5=0.64 
(5=0.36 


X 

ll 

-t] 

III 

ii 

^QNM jj 
. , . ■ i 


\ \ \ 

_ : v — * \- 1 

x \ \ 
s \ \ \ 

\ \ \ \ " 
\ \ \ \ 

\ \ \ N 



Mco 



FIG. 2: Energy spectrum for I = 2 and different values of j3. 
Horizontal lines are the corresponding ZFL-PP predictions, 
vertical lines are the QNM frequencies of the final BH. 



and the final ringdown tail. These seem to be universal 
properties of collisions involving BHs and were observed 
in the past in different settings [2(| The start of 

ringdown, roughly associated with the absolute maxima 
^20 ak m IV^oli occurs ~ 15M after the CAH formation, 
independently of 7. Except for a small neighborhood 
around 7 ~ 1, the maximal wave amplitude |^fo a I m " 
creases monotonically with the boost factor. The small 
dip in the wave amplitude for small, but non-zero veloc- 
ities has been seen before both in numerical simulations 
and analytic predictions (22J. For moderate boosts, we 
observe the absolute maxima in tp20 to be well approxi- 
mated by |Mr^ eak | « 0.26 + 0.48 7 ~ 2 [1/4 + log(l/2 7 )] 
[cf. Eq. ([3]) below] . The peak amplitude in the waveform 
h 2 o is roughly h^ k ~ ^20 7 W QNM' where wqnm is the 
lowest ringdown frequency for the mode [23j |. 

Figure [2] shows the energy spectrum JTJ) for collisions 
with different CM energy. For large CM energies, the 
spectrum is nearly flat up to some cutoff frequency. A flat 
spectrum is predicted by the ZFL and PP approaches, as 
indicated by the dotted lines in the figure. The cutoff fre- 
quency is well approximated by the least-damped QNM 
of the final hole, marked by a vertical line. The spectrum 
increases at small frequencies because of initial data con- 



tamination and finite-distance effects. 

Our numerical results indicate that the peak luminos- 
ity ([2]) is attained approximately 10M after the CAH 
formation. The peak luminosity is about 5 x 10~ 3 for 
/3 = 0.9, and may be as large as 10~ 2 as 7 — > 00. 
Restoring units, we get 10 _2 c 5 /G ~ 3.6 x 10 57 ergs~ 1 , 
the largest luminosity from a BH merger known to date. 
This is two orders of magnitude larger than for the infall 
from rest of two equal mass BHs, and one order of mag- 
nitude larger than for the inspiral of equal mass binaries. 
Nevertheless, it is still two orders of magnitude below the 
universal limit suggested by Dyson, dE/dt < 1 (24| . 

The total energy E radiated as a function of boost pa- 
rameter is shown in Fig. [3] Error bars on the radiated 
energies are determined as described in Sec. II. We have 
verified that E calculated from the radiation @ is con- 
sistent with alternative estimates obtained by directly 
measuring the mass of the final hole from the CAH prop- 
erties, and by using the ringdown frequency to estimate 
the mass of the final hole [23j . The ZFL predicts the fol- 
lowing functional form for the total radiated energy as a 
function of CM boost 7: 



E -w ( 1 ± 272 1 ^ - ^ log t ± V^ 1 ) 1 m 
M ~ y 2 7 2 + 2 7 3 V^T ) 

The quantity E^ is some unknown cutoff parameter, 
which is also the total fraction of energy radiated as 
7 — > 00. By fitting Eq. ^ to the numerical data we 
obtain = 0.14 ±0.03. The ZFL is a perturbative cal- 
culation about u> = 0, and its validity for our scattering 
problem is not obvious. However, given the good agree- 
ment with our numerical results in the kinetic-energy 
dominated regime 7 > 2, the extrapolation procedure 
should provide a reasonably accurate estimate for E^ . 

With regard to the multipolar contributions of the ra- 
diated energy, we find that E4 is at least one order of 
magnitude smaller than E2 for slow-motion collisions. 
This observation is consistent with the PP results for 
an infall from rest [2fj| , which predict an exponential de- 
crease of Ei with I. For larger boosts the ZFL and PP 
approach predict a strong increase in the relative contri- 
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TABLE I: Relative multipolar contribution (in %) and, in 
parentheses, the ZFL prediction. 

13 0.64 0.75 0.82 0.86 0.90 0.94 
E 4 /E 2 1.0(1.4) 2.4(3.4) 3.9(5.4) 5.0(7.3) 7.3(10) 11(14) 
WE e /E 2 0.2(0.3) 1.1(1.5) 2.1(4.0) 4.2(7.5) 11(16) 33(30) 



bution of higher multipoles, with Ei ~ M/l 2 as 7 — > 00. 
Our numerical results are in reasonable agreement with 
these calculations, as demonstrated in Table [U The dis- 
crepancies still present are due to the relatively large un- 
certainties in the energy carried by higher multipoles and 
to the breakdown of the ZFL prediction for small boosts. 
IV. Conclusions. In 1971, Hawking [25| placed an up- 
per limit of 29% on the total energy radiated when two 
BHs, initially at rest, coalesce. Numerical simulations 
of Einstein's equations [2l[ later showed that the true 
value is around 0.1% — two orders of magnitude smaller 
than Hawking's bound. Using a similar area theorem ar- 
gument, Penrose [3] derived an upper bound of 29% for 
ultra-relativistic head-on collisions (that the numerical 
values of the two bounds agree is apparently just a coin- 
cidence). Here we have presented results indicating that 
the answer in the high-energy limit is 0.14 ±0.03, slightly 
less than a factor of 2 of Penrose's bound, though quite 



close to the estimate of D'Eath and Payne computed us- 
ing perturbative techniques [26| . Even though our calcu- 
lations are in 4D, a consequence of this to searches for 
BH formation at the LHC is a warning that estimates of 
the "missing energy" based upon trapped surface calcu- 
lations could significantly overestimate this effect. 

This long overdue study represents an important step 
towards a full understanding of high-energy BH colli- 
sions. More accurate evolutions using significantly larger 
boosts are mainly inhibited by the junk radiation in the 
initial data. More work is also needed to study scat- 
tering with non-zero impact parameter, unequal masses 
and non-zero spins. For applications to LHC and RHIC 
physics, including the effects of extra dimensions, charge 
and AdS asymptotics (for RHIC) will be necessary. 
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